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Abstract 

The fluctuation properties of nuclear giant resonance spectra are studied 

in the presence of continuum decay. The subspace of quasi-bound states is 

specified by one-particle one-hole and two-particle two-hole excitations and 

the continuum coupling is generated by a scattering ensemble. It is found 

that, with increasing number of open channels, the real parts of the complex 

eigenvalues quickly decorrelate. This appears to be related to the transition 

from power-law to exponential time behavior of the survival probability of an 

initially non-stationary state. 
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Level fluctuations, measured in terms of the nearest-neighbor-spacing-distribution 
(NNSD) and the Aa-statistics, provide a commonly accepted tool for studying the quan- 
tum interplay between regular and chaotic dynamics. The standard treatment is restricted 
to bound states while, in many cases, the excited states are resonances embedded in the 
continuum. Already a generalization of the standard two-level repulsion theorem [|T| to res- 
onances shows that this may significantly modify the correlations between the states. 
Generically, chaotic dynamics leads to level repulsion but the presence of the continuum 
(open system), is expected P| to wash out the repulsion between the resonance energies. On 
the other hand, the lack of correlations between levels is normally interpreted as a manifes- 
tation of regular dynamics. It thus seems necessary to explore, on a fully quantitative level, 
what is the nature of the weakening of the repulsion due to openness and how it modifies 
the fluctuation characteristics. 

The most practical way for describing an irreversible decay into the continuum is based 
on a scattering ensemble of non-hermitian random matrices [Q. Such a treatment follows 
naturally from the projection-operator technique ^ in which the subspace of asymptotically 
decaying states is formally eliminated. The resulting non-hermitian Hamiltonian 

n = H-'-w (1) 

acts in the space of quasi-bound states and the coupling to the continuum is accounted for 
by the anti-hermitian operator W. Unitarity of the scattering matrix imposes on W the 
following factorization condition: 

W = AA^. (2) 

For an open quantum system with quasi-bound states, {i = 1, A^) which decay into 
k open channels a [a = 1, /c), the N x k matrix A = {^"} denotes the amplitudes for 
connecting the states \i) to the reaction channels a. The diagonalization of 7i in the basis \i) 
yields A^ quasi-stationary states with complex eigenvalues Sj = Ej —iTj/2, whose imaginary 
parts correspond to the 'escape width'. The factorization of W guarantees that > 0. An 
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interesting effect - due to the separable form of - is that, in the strong-couphng 
hmit (W ^ H), one observes a segregation of the states: k states accumulate most the total 
width, r = Z^jTj, while the remaining N — k states have nearly vanishing widths (they 
become 'enslaved' 0). 

For systems, such as the atomic nucleus, whose dynamics is expected to be classically 
chaotic, it is natural to consider the hermitian- and the anti-hermitian parts of Ti. to be 
statistically independent 0. Furthermore, the real and symmetric N x N matrix H can be 
modeled @] as a member of the Gaussian orthogonal ensemble (GOE) of random matrices 
|0J^. For large the matrix elements of H obey the following pair contraction formula: 

{Hii'Hjj') = -^{SijSiij' + 5ij'5i'j) (3) 

in the sense of GOE averaging. The constant a is related to the mean level spacing, D = 
2a/ N. 

For a general Gaussian ensemble of complex random matrices 7i an analogous con- 
traction formula for {Tiijl-iii ji) is obtained which implies that the real and imaginary parts 
of Ti commute on average. Consequently, the two hypersurfaces, representing the real and 
imaginary parts of the energy lie in orthogonal subspaces []TU[. This, for sufficiently large N , 



may produce decorrelated spectra as seen from either the real or imaginary axes, in spite of 
a cubic repulsion on the complex plane. 

However, this general Gaussian ensemble of complex random matrices is not applicable 
in the present case because of S'-matrix unitarity. Instead, the anti-hermitian part of Ti 
is determined by the amplitudes A"^ via Eq. (^. Based on the GOE character of internal 
dynamics and orthogonal invariance arguments |^ the amplitudes A" can be assumed to be 
Gaussian distributed. The corresponding correlator reads: 

(^•4) = 1^1"^"%^ = (4) 

implying that the average trace is (TrW) = Sa7". The diagonal elements Wu = S^^^(A")^ 
are then positive, statistically independent and obey a Xfc-square distribution. 



Unlike the amplitudes the matrix elements of W are not statistically independent, 
however. The number of independent random parameters, Nk — ^k{k — 1) for k < N, is 
reduced by the second term as a consequence of the rotational invariance of Wij = S^^^^A^A" 
(the scalar product between N /c-dimensional vectors Aj in the channel space). Only for 
k = N the correlations in W are specified by ^N{N — 1) parameters, as for the GOE. Thus 
a decorrelation of the projected spectra may result. In most realistic cases, however, the 
number of open channels k is smaller than A^. To assess the dependence on the number 
of open channels we perform a systematic numerical study of the spectral correlations as a 
function of k. 

Since the nuclear interaction is predominantly two body in nature, the matrix represen- 
tation of the nuclear Hamiltonian should be related to the so-called 'embedded' Gaussian 
orthogonal ensemble (EGOE) rather than the GOE. Therefore, to make our study real- 
istic from the nuclear physics point of view, we generate the hermitian part of Ti from the 



model in ref. |Tl| instead of using a GOE random ensemble. The Hamiltonian includes a 
mean-field part and a zero-range and density-dependent two-body interaction. The matrix 
representation of H is expressed in the basis of one-particle one-hole (Iplh) and two-particle 
two-hole (2p2h) excitations generated by the mean-field part and by discretizing the con- 
tinuum [0]. The spectral fluctuations of the corresponding real eigenvalues, measured in 



terms of the NNSD and A3, coincide with those of the GOE |]1T|, even though significant 
deviations from the Gaussian distribution of the matrix elements are found 

Because of time-reversal invariance the anti-hermitian part of Ti is generated by a Gaus- 
sian ensemble of real amplitudes A'} with correlator (^, where 7° = 1, i.e. we assume 
that all channels are equivalent and the strength of the external coupling is comparable to 
the internal one. In the specific calculations presented below, we select quadrupole exci- 
tations (J'^ = 2"*") in ^°Ca. To ensure acceptable statistics, in the quasi-bound-state space 
all Iplh and 2p2h states up to an excitation energy of 40 MeV are included. This yields a 
1661 X 1661 Hamiltonian matrix. Fig. 1 shows the resulting eigenvalue distribution on the 
complex energy plane for an increasing number k of open channels. For A; = 10 the majority 



of the energies lie very close to the real axis and only a few states acquire a significant 
width which is a trace of the 'collective synchronization' discussed in ref. Increasing 
k, the distribution becomes more uniform and the width of the empty strip between the 
cloud of eigenvalues and the real axis widens. This is understandable as is equal to the 
'correlation width' which describes the asymptotic behavior of the decay process . 



The NNSD on the plane can be determined by calculating the normalized distances 
Si = diPniSiY^"^, where di stands for the Euclidean distance between the eigenvalue Si and 
its nearest neighbor, and Pn{Si) for the local density of eigenvalues determined from n 
nearest neighbors of £i. Similarly as in ref. [0, the choice = 10 turns out satisfactory 



and guarantees stability. The numerical results are compared to the Poisson distribution 
P{s) = (7r/2)s exp(— 7rs^/4) (dashed lines in the rh column of Fig. 1), which shows linear 
repulsion on the plane, and to the P{s) = (8l7r^/128)s'^ exp(— 97rs^/16) with cubic repulsion 
(solid lines). The latter gives a good description for the NNSD of symmetric Gaussian 



random matrices [|T2| |[T^ and, for a large number of open channels, also fit our numerical 
results nicely. For a few open channels (upper right part of Fig 1.) we see a weaker then 
cubic repulsion, however . 

Now we come to the central point namely the fluctuation properties of the real parts Ei 
of the energy eigenvalues. The corresponding NNSD and As-statistics are shown in Fig. 2. It 
is well known that, without coupling the continuum, the spectra show GOE characteristics 



for both measures |jTT]]. However, for many open channels a decorrelation takes place. In 
fact, for large k the results are well reproduced by a Poissonian shape of the NNSD (lower 
left part of Fig. 2). Quite surprisingly, this even holds for k/N oi a. few percent (middle left 
part of Fig. 2). Already for ten open channels {k/N = 6 * 10^^), there is a visible deviation 
from the Wigner distribution (upper left part of Fig. 2). These numerical observations lead 
to the conclusion that the appropriate way of describing these deviations is to consider 
superpositions of Wigner and Poisson distributions rather than Wigner and Gaussian . 

The longer-range correlations (spectral rigidity) expressed by the As-statistics show a 
similar tendency, although the transition is somewhat slower. In addition, as is seen in 



Fig. 2, the transition region Lmax from GOE to Poissonian characteristics is restricted to 
about 10 normahzed distance units. This appears to be consistent with the findings in 
p!8| for hermitian separable problems, where L^ax increases with increasing length of the 
string of eigenvalues. In the present case the string is comparatively short. On a more 
formal level the As-statistics is known to be non- universal above a certain L^ax- For 
systems with a known classical limit, Lmax is determined by the inverse of the period of 
the shortest periodic orbits. We wish to mention, without showing the results explicitly, 
that an analogous analysis for the imaginary parts of £i show Poissonian fiuctuations for 
any number of the open channels. This asymmetry in the statistical properties of Ej and 
Tj is related to the different properties of the real and imaginary parts of ?i, especially for 
smaller values of k. 

Another way of understanding the decorrelation of the resonance energies due to the 
presence of continuum decay comes from the relation between the wave-packet dynamics 
and the stationary states [^. The latter can be obtained via the Fourier transform of the 
time evolution of a generic wave packet. For a bound-state problem such a wave packet 
resides in the interaction region forever and thus, the structure of the corresponding phase 
space can be resolved with arbitrary accuracy. Consequently, for a chaotic system, the whole 
complexity (delocalization, random nodal pattern, scars, etc.) of stationary states can be 
reproduced. Coupling to the continuum, sets a limit for this process, however. As time 
progresses, the wave packet will leak out of the interaction region and makes it impossible 
to resolve all details of the dynamics. As a result the wave functions, projected onto the 
interaction region, look more regular than their counterparts in a closed system. The leakage 
is expected to occur faster with increasing k. A quantititive measure of the speed is the 
survival probability P{t) of a randomly chosen wave packet initially localized in the 
interaction region. As a convenient and experimentally motivated choice we consider a state 
excited by the isovector quadrupole operator = -F|0)). When expanded \F) involves all 
the eigenstates \xi) of Ti and 
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N 

Pit) = |(F(0)|F(t))r = I E(0|^lx.)(x,l^|0)e^^^*/^P (5) 

(for a complex symmetric matrix the left and right eigenvectors are the same). In the absence 
of continuum coupling, P{t) remains constant (on average) after a rapid initial dephasing 
due to the non-stationarity of \F) [^. For an open system, on the other hand, a decay 
of P(t) is to be expected. The most interesting feature is the dependence of the decay law 
on the number of open channels: For a small k the decay is very slow and well represented 
by a power-law {P{t) ~ t^^)- For k = 1 we find z ~ —1/2, in reasonable agreement 
with the estimates of ref. |2^. As k increases z grows very fast and, for k > 100, P{t) 



drops exponentially on long time scales, i.e. P{t) ~ exp(— ?7t), with the decay constant r) 
growing rapidly with k (Fig. 3). These observations go in parallel with the classical picture 
of open phase space phenomena such as a chaotic scattering [23|: For a small number of 



the open channels the decay is governed by a power-law. This is associated with larger 
fractal dimensions of the set of singularities generating chaotic behavior than for many open 
channel cases which lead to an exponential decay. 

In summary, the numerical analysis presented in this work shows that GOE correlated 
spectra of quasi-bound states become fully decorrelated in the presence of continuum cou- 
pling and when the number of open channels is large. This transition is accompanied by a 
change of the decay properties of the average survival probability of a non-stationary wave 
packet, turning from power-law to exponential. This appears to be consistent with the 
semiclassical relation between the time-dependence of P{t) and the structure of the res- 
onances. An exponential behavior of P(t) corresponds to the region of strongly overlapping 
resonances (Ericson fluctuations [^), while the power-law decay, with small power indices 



z |22l , corresponds to isolated resonances, and it is this isolation which preserves the original 
fluctuations. 
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Figure Captions 

Figure 1: Left column: The eigenvalue distribution of the non-hermitian Hamiltonian Ti defined 
in Eq. (1) for different number k of open channels. The hermitian part H is chosen 
as the Hamiltonian of [12] while the anti-hermitian part W is given by Eq. (2) tak- 
ing the amplitudes A as members of the Gaussian ensemble [4]. Right column: the 
corresponding NNSD on the complex plane. 

Figure 2: The NNSD [Ihs) and the A3 statistics (ris) of the real parts Ei for energy eigenvalues 
of Ti. and different number k of open channels. 

Figure 3: The time dependence of the survival probability P{t) of a wave packet, initialized by 
the isovector quadrupole operator, for various numbers of open channels. 
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